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Abstract: The so-called spectral dimension is a scale-dependent number associated 
with both geometries and field theories that has recently attracted much attention, 
driven largely, though not exclusively, by investigations of causal dynamical triangu- 
lations (CDT) and Hofava gravity as possible candidates for quantum gravity. We 
advocate the use of the spectral dimension as a probe for the kinematics of these (and 
other) systems in the region where spacetime curvature is small, and the manifold is flat 
to a good approximation. In particular, we show how to assign a spectral dimension (as 
a function of so-called diffusion time) to any arbitrarily specified dispersion relation. 
We also analyze the fundamental properties of spectral dimension using extensions of 
the usual Seeley-DeWitt and Feynman expansions, and by saddle point techniques. 
The spectral dimension turns out to be a useful, robust and powerful probe, not only 
of geometry, but also of kinematics. 
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1 Introduction 

The spectral dimension has been proposed as possible observable characterizing the 
geometry in discrete quantum gravity approaches [1]. It can be viewed as an effective 
notion of dimension defined though a fictitious diffusion process on the discrete geom- 
etry. In simple terms (see reference [1] for more details), the diffusion process can be 
thought as a stochastic random walk and the spectral dimension is defined in terms of 
the average return propability P(s) of the fictitious diffusion as 



with s being the (fictitious) diffusion time. As opposed to the topological dimension d, 
the spectral dimension ds(s) need not be an integer and furthermore is scale-dependent. 

The concept of spectral dimension has attracted considerable interest in causal 
dynamical triangulations (CDT) [1-6], since finding meaningful observables in discrete 
theories is not a trivial task. (See also references [7-9] where the spectral dimension is 
used in noncommutative spacetimes, Euclidean dynamical triangulations and random 
combs respectively.) Such observables can provide the much needed connection with 
continuum theories (at suitable limits). Indeed the spectral dimension can be defined 
in the continuum quantum gravity models as well and be used to characterize and 
understand their short-distance behaviour (see references [10-14]). Furthermore, it was 
shown in reference [10] that both CDT and Hofava gravity seem to lead to a spectral 
dimension of 2 in the ultraviolet limit, while the value of the spectral dimension matches 
that of the topological dimension at larger length scales. These first results suggest that 
the spectral dimension can be a useful potential link between discrete and continuum 
theories. 

Clearly if one is to utilize the spectral dimension in this way, it is crucial to first 
get a deep understanding one what kind of information it actually carries in continuum 
theories and also develop a set of tools for extracting this information. Our goal here 
is to make a first step in this direction. 

Given the definition of a spectral dimension, it is actually quite intuitive (and will 
become clearer shortly) that, in continuum theories, its behaviour is scale-dependent 
and one can expect that: 
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1. At long distances, scales which are comparable to, or larger than, the radius of 
curvature, the spectral dimension will depend on the details of the geometry and 
the presence gravitational sources. At these scales the deviation of the spectral 
dimension from the topological dimension can be essentially attributed to curva- 
ture effects. Alternatively, at these scales the spectral dimension probes aspects 
of the geometry. 

2. At intermediate scales, large compared to the Planck scale but small compared 
to the radius of curvature, spacetime is effectively flat, and one should obtain a 
spectral dimension whose value matches that of the usual topological dimension 
of spacetime. 

3. At short distances, (typically Planckian distances, but sometimes one is instead 
interested in the near-horizon limit of a black hole), curvature effects can certainly 
be neglected, and the deviation of the spectral dimension from the topological 
dimension should be attributed to different effects. 

We will not be directly interested in regime 1 for the purposes of the current analysis. 
That is, we will not attempt any thorough analysis of how the spectral dimension can 
be used as a probe of geometry or how curvature affects the spectral dimension (though 
we will comment upon it briefly). Instead, in this article we shall focus on the flow 
of the spectral dimension from flat spacetime, which is for our purposes the infrared 
limit, into the ultraviolet limit (regime 2 above to regime 3 above). The corresponding 
length scales are of particular interest in quantum gravity candidates. We shall develop 
a number of theoretical results, and complement them with analytic and computational 
tools that allow for both exact and approximate evaluation of the spectral dimension 
in a number of interesting situations. 

In particular, we shall demonstrate that the spectral dimension is not necessarily 
intrinsically geometric. At scales small enough for curvature effects to be negligible 
its deviations from the topological dimension it actually becomes an analytic property 
of the differential operator (ultimately equivalent to a dispersion relation) that one 
is using as input to define the generalized diffusion process (as already suggested by 
the analysis of reference [10] for the particular case of Hofava gravity). In turn, this 
operator acts as the propagator of some dynamical degree of freedom in flat space. In 
this sense, the spectral dimension acts, at suitable scales, as a probe of the kinematics 
of the particular degree of freedom. Within this framework we will show that, given a 
specified topological dimension d, it is possible to define a scale-dependent notion of a 
spectral dimension ds(s) for any arbitrary dispersion relation u = Q(k). We shall also 
derive a formal inversion process, so that from the spectral dimension, (as a function 
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of s), and the topological dimension d, one can in principle reconstruct the dispersion 
relation. 

For some specific dispersion relations the situation is sufficiently simple that it is 
possible to explicitly evaluate the relevant integrals, and so explicitly determine the 
spectral dimension in terms of standard special functions (such as Bessel functions and 
[Gaussian] error functions). More generally, we shall develop a number of analytic tools, 
(using a generalization of the Seeley-DeWitt expansion in the ultraviolet, and a gener- 
alization of the Feynman expansion in the infrared), that give us good analytic control 
of these two limits. We also develop an approximation based on saddle point tech- 
niques to express the spectral dimension in terms of the phase and group velocities of 
the underlying dispersion relation — this approximation agrees with asymptotic results 
in the ultraviolet and infrared, and can also be used in the intermediate regime. (When 
saddle point and exact integral techniques overlap, the saddle point approximation will 
be seen to be a tolerably good global fit to the exact result.) 

Apart from those issues of relevance to the quantum gravity and field theory com- 
munities, we additionally point out that the formalism derived in this article may also 
be of relevance in completely non-relativistic situations (such as the propagation of 
surface waves on liquid-gas interfaces.) Overall, we wish to advocate the use of spectral 
dimension as a general, powerful, and robust probe of the kinematics (specifically the 
dispersion relation) of whatever system one has under consideration. 

2 Theoretical framework 
2.1 Dispersion relation 

Suppose we are working in flat (d + 1) dimensional spacetime and have a (in general 
Lorentz-violating) dispersion relation 

u = n(k). (2.i) 

Now define 

f{k 2 ) = Q(k) 2 . (2.2) 

This dispersion relation can always be viewed as being completely specified by the 
solutions of the (in general Lorentz-violating) differential equation 

D L ^ = {-«9 2 -/(-V 2 )}$ = 0. (2.3) 

The reason that we choose the differential operator to be second order in physical time 
t is a purely pragmatic one — with such a choice the differential equation encoding the 
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dispersion relation can typically be derived from a ghost free non-interacting Lagrangian 

jgf = I = I $ - /(-V 2 )} $. (2.4) 

2.2 Spectral dimension 

First we Wick-rotate in physical time t to consider the Euclideanized differential oper- 
ator De in (d + 1) topological dimensions 

= {-<9 2 + /(-V 2 )}$. (2.5) 

Now consider the generalized diffusion process 

^ + D E ^p(x,x',s) = 0; p(x,x',0) = 5 d+1 (x-x'), (2.6) 

where D E is the differential operator defined above. We set x = (t,x), and s is best 
thought of as an auxiliary diffusion time. (That is, s is a "fake time" , a Lorentz- violating 
generalization of Schwinger-De Witt proper time — as in the usual Schwinger-De Witt 
proper-time regularization of spacetime propagators.) We shall soon see that s can also 
be thought of as the scale at which one is probing the diffusion process — we shall see 
that s — y probes the ultraviolet, while s — > oo probes the infrared. The general 
solution of equation (2.6) is given by 

/c\ d k c\u) -> 
exp{i[* ■ (x - f ) + u(t - t')}} exp{-s[co 2 + f(k 2 )]}, (2.7) 

as can easily be verified by differentiation and by noting that the boundary condition 
is satisfied at s — 0. The quantity p(x, x' , s) is often called the heat kernel and can 
be thought of as the probability density that a particle originally at x' at s = will 
diffuse to x in diffusion time s. Now consider the average return probability density 
P(s) (also called the diagonal part of the heat kernel): 

/c\ d k c\h) 
^-^ex P {-s[u 2 + f(k 2 )}}, (2.8) 

(We are working in flat space where, due to translational invariance the average return 
probability equals the return probability.) The spectral dimension is defined as 

d s (s) = -2 dl " P(s) = -2^#M. (2.9) 
d In s a In s 

The reason for this definition is that if we consider the particularly simple differential 
operator D — □ = —V 2 , then in (d + 1) topological dimensions then one can show that 

P_ y2 ( s ) =p_ v2 (x,x;s) oc s - (d+1)/2 (2.10) 
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so that in this specific case 



d 



S;-VHSJ 



> dlnP - v2(g) = _ 2 dlnp-y 2 (*,x; g ) =d+1 
d In s d In s 



That is: The spectral dimension d s -D E ( s ) °f the complicated diffusion process governed 
by the differential operator De is the (fictitious) topological dimension of that sim- 
ple diffusion process governed by —V 2 that most closely approximates Pd e (s) at the 
indicated value of the diffusion time s. We emphasise that this is a definition, not a 
theorem. This definition is useful only insofar as it leads to interesting results. 



2.3 From dispersion relation to spectral dimension 

Factorize P (s) into "time" and "space" contributions 

d d k 



P s) 



(2.12) 



The first integral is elementary, so one has 



P(s) 



.V4 



7TS 



X 



d d k 



(27T) 



T exp{- S /(^)} 



(2.13) 



Dropping irrelevant constants of proportionality (they would in any case drop out in 
the calculation of the spectral dimension) we see 



P(s) oc 4= / ^- 1 exp{-s/(P)}dA;. 
V s J 



(2.14) 



That is 



lnP(s) = -^lns + ln J k d ~ l exp{-s f(k 2 )}dk + C. (2.15) 

(Note that C will be used to denote a generic constant, when we re- use the symbol 
below it will not necessarily represent the same constant.) This expression lets you 
calculate the exact spectral dimension as 



d s (s) = 1 + 2s 



J f(k 2 ) k*- 1 exp{-s/(£; 2 )}d£; 
JW 1 exp{-s/(fc 2 ) }dk ' 



(2.16) 



Note that the leading "1" comes from the fact that the physical time part of the 
differential operator De is trivial — it is just d 2 . All the complications come from the 
d spatial dimensions. If we now write the dispersion relation as 

oo 2 = f(k 2 ) = n(k) 2 , (2.17) 
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then, still as an exact result: 



d s (s) = l + 2s 



jfl(k) 2 k 11 - 1 exp{-sQ{k) 2 }dk 
JW 1 exp{-sQ{k) 2 } dk 



(2.18) 



It is now useful to define a partition function 



Z(s) = / k d ~ l exp{-sf(k 2 )}dk = / k d ~ l exp{-stt(k) 2 }dk. (2.19) 



Now by construction Z(s) oc P{s), so on the one hand it carries no extra information. 
On the other hand it lets us formally rewrite the spectral dimension as an ensemble 
average: 



Note that this partition function encodes relatively simple information concerning the 
dispersion relation of the specific degree of freedom under consideration. It is by no 
means the partition function of the entire system. 

This completes the first (and simplest) part of the programme: From any arbitrary 
dispersion relation u = Q(k), and specified topological dimension d, we have seen how 
to construct suitable differential operators Di amd De that encode this dispersion 
relation, and how to use the generalized diffusion process associated with the differential 
operator De in order to define the corresponding spectral dimension ds(s). 

2.4 From spectral dimension to dispersion relation 

The definitions given above allow one to determine the spectral dimension for a given 
dispersion relation. This is the usual case when one deals with a continuum theory and 
can linearize the field equations in order to determine the dispersion relation. In such 
cases the spectral dimension acts as an additional probe of the kinematics described 
by this dispersion relation (and potentially the geometry, see also section 2.5.2 below). 

There are, however, cases where one can calculate the spectral dimension by other 
techniques and it then becomes useful to determine a corresponding (effective) disper- 
sion relation. A typical example arises in discrete theories, such as CDT, where one 
can define and calculate the spectral dimension directly from a discretized diffusion 
process [1, 6]. Clearly, in such theories the spectral dimension is not just an alternative 
to studying directly the propagator. Since the latter is simply not available, the former 
becomes a far more essential tool for obtaining information about the theory. Deter- 
mining a corresponding (effective) dispersion relation could facilitate the comparison 
to some candidate continuum effective theory, see reference [15] for an example. 




(2.20) 
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2.4.1 Dispersion relation as inverse Laplace transform 

Formally one can invert the process above and from the spectral dimension (as a func- 
tion of s), plus the specified topological dimension, reconstruct the dispersion relation. 
This is best done in two steps. From 

dlnZ(s) d s (s)-l 



ds 2s 

we have 



(2.21) 



_1 p d 8 ( 8 ) l d _ 



■S|) 



That is, up so some constant of proportionality Z(sq), we can certainly reconstruct the 
partition function Z(s) from the spectral dimension. Can we take this any further? 
Can we reconstruct the dispersion relation Q(k) from the spectral dimension ds{s)7 
After some mathematical manipulations one can bring equation (2.19) in the form 

ZW = ^~exp { -^}^»d M . (2.23) 
Performing an integration by parts yields 

POO J 

/ expi-sfi 2 } {k{Q) d } d[Q 2 } = -Z(s). (2.24) 
Jo s 

This has the form of a Laplace transform, in the variable Q 2 , of the function k(Q) d . 
Implementing the inverse Laplace transform via a complex integration we have 

k(n) d = — / - Z(s) e n2s ds, (2.25) 
2m J c s 

where C is an appropriate contour in the complex plane and Z(s) is given by equa- 
tion (2.22). Therefore, in principle one can calculate the dispersion relation when the 
spectral dimension is analytically known on the complex plane as a function of s. 

2.4.2 Practical inversion techniques 

Though it settles an important issue of principle, the method of calculating the disper- 
sion relation from the spectral dimension presented above is usually of little practical 
value. This is because, in either discrete or continuum theories, when one might be in- 
terested in determining a corresponding dispersion relation from a spectral dimension, 
the latter is typically not known analytically as a function of s. Hence, an algorithm 
that requires analytic continuation of the spectral dimension to complex values of s is 
impractible. This problem can be circumvented in various ways: 
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• Direct inversion via semi-analytic model building: If one has somehow, (for in- 
stance, from some numerical simulation followed by some data-fitting), built a 
semi-analytic model for the spectral dimension ds{s), then integrating to find 
the partition function Z(s) is typically easy. This semi-analytic model partition 
function Z(s) can be analytically continued to the complex plane, and sometimes 
one is lucky — one may encounter a model for the partition function for which 
the inverse Laplace transform is explicitly known. Even when this works, we 
emphasise that one is analytically continuing the model, not the data. 

• Direct inversion via infinite differentiation: There is a little-known method due 
to Post [16], see also Bryan [17], that allows for inversion of Laplace transforms by 
taking arbitrarily high derivatives. Specifically, if G(s) is the Laplace transform 
of g(t) then 



This algorithm may not always be practical, since one needs arbitrarily high 
derivatives. Even if not always practical, it again settles an important issue of 
principle - - knowledge of the spectral dimension ds{s) in principle allows one 
to reconstruct an equivalent dispersion relation Q(k). 

• Indirect inversion via nonlinear regression: In contrast, a practical algorithm for 
extracting (or rather, approximating) the dispersion relation is to build some 
theoretically appealing model for Q(k) with several adjustable parameters, and 
to then use the technique of nonlinear regression to find a best fit within that 
class of dispersion relations. In a companion paper [15] we have applied this idea 
to the measured spectral dimension function arising in CDT models of quantum 
gravity. 

2.5 Some generalizations 

Let us now indicate some possible generalizations and extensions of this formalism. 
2.5.1 Beyond second order in physical time derivatives 

If one is (perhaps foolishly) willing to play with ghost fields then there is no need to 
restrict to second order derivatives in physical time. Consider the operator 




(2.26) 




(2.27) 



The relevant dispersion relation is now only implicitly defined by 




(2.28) 
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On the other hand the heat kernel becomes 

id 



exp{i[k ■ (x - x') + u(t - t')}} exp{-s[f(u 2 , k 2 )}}, (2.29) 
and the diagonal part of the heat kernel specializes to 

/r\ d k r\t,! 
^-^exp{- S /(u, 2 ,A; 2 )}. (2.30) 

So the spectral dimension is 

, M „ / k 2 ) exp{- g /(a; 2 , k 2 )} da; fc^ 1 dfc 0/ff3j *„ f0 ^ 

= 25 Jexp{- g/ (c 2 ,P)}da;^ dfe = 2s(/( " ' * ))s (2 " 31) 

and depends on the whole of the function f(u 2 ,k 2 ), while the dispersion relation is 
given by the solution of equation (2.28) and uses only limited information about the 
function f(cu 2 , k 2 ). 

If the differential operator can be put in the form 

^$^{/ 1 (-9 2 ) + / 2 (-V 2 )}<l>, (2.32) 

then the space and time contributions to the spectral dimension become separable and 
one has 

, n _ 9 M(^)exp{-.A(a, 2 )}da; J/ 2 (fc 2 ) fc^ 1 exp{-s / 2 (fc 2 )} dfc 
51 j " /exp^aAM}^ J^-i exp{- S / 2 (fc 2 )}dfc " 1 66} 

Clearly, when fi(—d 2 ) = —d 2 , one recovers the result given in equation (2.16). 
2.5.2 Curved spacetime 

A perhaps more interesting question is to ask what happens if spacetime is allowed 
to be curved, in addition to having a nontrivial dispersion relation. Extending the 
definition of the spectral dimension to dispersion relations describing propagators in a 
given curved background is actually rather straightforward: De in equation (2.6) will 
be a differential operator in curved space, but the solution of this equation can still be 
used to calculate the return probability. The only subtlety is that when one has to deal 
with spaces which do not respect translational invariance, the return probability would 
be space dependent, P(x, s). Thus, so long as one does not want the spectral dimension 
to depend on the spacetime point, one has to actually use the average return probability 
in the definition of equation (2.9) (as was mentioned already in the the introduction). 
An Arnowitt-Desser-Misner decomposition provides a natural setting for constructing 
differential operators which include curvature and at the same time lead to nontrivial 
dispersion relations even around maximally symmetric spaces. 
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2.6 Space versus spacetime 



The ds{s) we have so far been considering is the spectral dimension of spacetime. If 
one is only interested in the spectral dimension of space, (as is usual for most classical 
[non-relativistic] physicists), then the appropriate thing to do is to simply drop the t 
coordinate and define 

Apace = -/(-V 2 ). (2.34) 
One would then consider the diffusion process 
d 



^ - /(-V 2 ) | P {x, x', s) = 0; p{x, x\ 0) = 5 d (x - x'). (2.35) 

In this case (now without the leading "1") we have 

ff(k 2 ) k d ~ l exp{-s f(k 2 )}dk 
ds — {s) = 23 Jk^ex P {-sf(k 2 )}dk ' (2 - 36) 



or equivalent ly 



/ il(k) 2 k d ~ x exp{-sn(k) 2 }dk 



ds ^ = 2s j^«p{-W}d* = 28 (m ) - ■ (2 ' 37) 

Note that (assuming physical time only shows up in the simple form d 2 ) the only 
difference between the spectral dimension of spacetime, and that of space itself, is the 
leading "1" . Furthermore, since t has been eliminated from the formalism, the diffusion 
time s can in this situation be identified with physical Newtonian time. 



3 Calculational techniques 

We shall now explore some calculational techniques for determining the spectral di- 
mension ds{s). In simple situations we can obtain exact analytic results, but more 
typically one has to resort to some form of approximation. We shall show that making 
use of asymptotic expansions can be fruitful in both the infrared and ultraviolet, and 
that saddle point techniques can also be employed — typically over the whole range of 
interest. 

In the following discussion we will also present several examples. For simplicity we 
will largely focus on polynomial dispersion relations, where 

uj 2 = Cw k 2w + c w+1 k 2w+2 + ■■■ + c z _ x k 2z ~ 2 + c z k 2z . (3.1) 
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Here w < z are integers potentially unrelated to d. As additional examples (better dealt 
with using asymptotic or saddle point techniques) we will also use rational polynomial 
dispersion relations of the form 

u - w~r ( } 

where p(k 2 ) and q(k 2 ) are polynomials in k 2 , and surface-wave dispersion relations that 
sometimes include trigonometric functions. 

The choice of these specific types of dispersion relations as examples is motivated 
more by mathematical simplicity and their demonstrative power, than by their origin 
or genericity. Nevertheless, all of them are common (to different extents). Polynomial 
dispersion relations are typical in (truncated) effective field theories that exhibit Lorentz 
violations. They also appear in projectable Hofava gravity [18-21]. (Projectability 
requires that the lapse function is a function of time only, N = N(t).) In the latter 
case, the dispersion relation is actually polynomial in k 2 with the highest power being 
(k 2 ) d , that is, 

uj 2 = Cl k 2 + c 2 k A + ■■■ + c d k 2d . (3.3) 

and with particularly pleasant renormalizability properties being associated with the 
choice z = d [18, 22, 23]. Rational polynomial dispersion relations are generic in more 
general, non-projectable Horava gravity [24, 25]. (See reference [26] for a brief review 
on Hofava gravity.) Finally, gravity-driven surface waves have dispersion relations 
that include trigonometric functions. In any case, the techniques we present here are 
applicable to more general dispersion relations. 

3.1 Exact results 

Useful exact results seem to be limited to special cases of the two-term dispersion 
relation 

u 2 = c w k 2w + c z k 2z . (3.4) 

More complicated dispersion relations, or even the general case above, typically require 
one to invoke some sort of series expansion or approximation technique. 

3.1.1 Bogoliubov dispersion relation 

The so-called Bogoliubov dispersion relation can conveniently be put in the form 



k 4 / k 2 

/(* 2 ) = * a + 2T^; n(*) = *Vi + 2r^- (3-5) 

(The factor of 2 d is there just to minimize irrelevant numerical factors in the results 
below.) 
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Two space dimensions: The partition function for the Bogoliubov dispersion rela- 
tion can be explicitly evaluated in terms of the error function 



Z{s) 



kexp 
K Fk 



sk 

exp (s-fT 2 ) 



1 + 



k' 2 



dk 



2 V s 

The spectral dimension is then explicitly evaluable as 



AK 2 J 
1 - erf f VsK 2 



(3.6) 
(3.7) 



d s (s) = 2- 2sK 2 + 2 



sK 2 



TV 



sK 2 



eriiVsK 2 ) 



Note the limits 



d s (s -»■ 0) = 2; d s {s -»• oo) = 3. 



(3.8) 



(3.9) 



Three space dimensions: In this situation the partition function for the Bogoliubov 
dispersion relation can be explicitly evaluated in terms of Bessel functions: 



zoo 



fc 2 exp <^ -sk 2 1 



k' 2 



8K'- 



dk 



= V2K 3 exp (sK 2 ) [K 3/4 (sK 2 ) - K 1/4 (sK 2 )] 
The spectral dimension is then explicitly evaluable as 



Note the limits 



dots) = 2- 4sK 2 - 



d s (s 0) 



K 3/A (sK 2 ) + K l/4 (sK 2 
2 [K 3/4 (sK 2 )-K 1/4 (sK 2 )\ 

5 

-; da(s — > oo) = 4. 



(3.10) 
(3.11) 

(3.12) 
(3.13) 



3.1.2 Quadratic plus sextic dispersion relation 

Consider the specific dispersion relation 



f(k 2 



1 + 



27K 4 



k f 



21K V 



Q(k) = k\ 1 



k 4 



27 K 4 ' 



(3.14) 



(The "27" is there just to minimize irrelevant numerical factors in the results below.) 
This dispersion relation is potentially interesting as can serve as a good approximation 
for both the high-momentum and low-momentum limit of Hof ava gravity [10] (it inter- 
polates between low-momentum Lorentz invariance and a z = 3 Lifshitz point at large 
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momentum). It is not, however, the exact dispersion relation of any of the versions of 
the model in 3+1 dimensions, (there is no k 4 term, and it is not a rational polynomial). 
The partition function can be explicitly evaluated in terms of Bessel functions 



Z(s) 



Arexp <^ -sAr 1 + 



k 4 



27K*) 



dk 



(3.15) 



= -t^K^s [J 1/3 (sK 2 ) J_ 2/3 (sK 2 ) + Y 1/3 (sK 2 ) Y. 2/3 (sK 2 )] . 
The spectral dimension is then explicitly evaluable as 



d s ( s ) = 2 + 2sK 2 
Note the limits 



J 1/3 (sK 2 ) 2 - J_ 2/3 (sK 2 ) 2 + y 1/3 (sK 2 ) 2 - F_ 2/3 (sK 
J 1/3 (sK 2 ) J_ 2/3 (sK 2 ) + Y 1/3 (sK 2 ) Y_ 2/3 (sK*) 



d s (s -»■ 0) = 2; ds(s -»■ oo) 



4. 



(3.16) 



(3.17) 



3.1.3 Quartic plus sextic dispersion relation 

As a final exact example, albeit one where the results are too messy to actually explicitly 
write down, in (3+1) dimensions consider the dispersion relation 

2k 6 



f(k 2 



3k* + 



K 2 



(3.18) 



The numerical factors are again chosen to keep expressions as simple as possible. The 
partition function can then be evaluated as a sum of generalized hyper-geometric func- 
tions, specifically as a sum of three 2 -F 2 (— sK 4 ) hyper-geometric functions. We have 



Z(s) 



2k® \ 



k z exp{-s {3k ;* + —j<jdk 

2-p2([l/4,3/4], [1/3,2/3]) (-S-^ 4 ) 



1 2nK 2 

12V^ 



(2s K 4 ) 1 / 3 2 F, 



2([7/12,13/12], [2/3,4/3]) \ 



-sK 



isr(j) 

16^ 



(2sK 4 ) 2 / 3 2 -F 2 ([n/i 2i i7/i2] 



[4/3,5/3]) ( 



-sK 4 



(3.19) 



The resulting spectral dimension is too messy to be worth writing down explicitly, but 
the limits are straightforward to evaluate: 

5 
2' 



d s (s -»■ 0) = 2; d s (s -> oo) 



(3.20) 
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3.1.4 Summary of exact results 

It should be obvious by now that even relatively simple dispersion relations, such as 

cu 2 = c w k 2w + c z k 2z , (3.21) 

lead to rather complicated expression for the spectral dimension. Therefore, it is inter- 
esting to develop techniques that could allow us to determine the spectral dimension 
to some required accuracy based on a suitable approximation. This will be the aim of 
the in the following sections. 

Apart from demonstrating the difficulty of determining the spectral dimension for 
a given dispersion relation, the exact results presented above provide also some insight 
into how the spectral dimension changes with the length scale. Given that s is a 
(fictitious) diffusion time, small values of s probe small length scales and large values 
of s probe large scales. This is the reason why we paid particular attention to the limits 
s — y and s — > oo. For dispersion relations of the form given above and whenever we 
have been able to explicitly carry out the exact calculation we have found that 

d s ( s -> 0) = 1 + -; d s (s -> oo) = 1 + — , (3.22) 

z w 

in (d + 1) topological dimensions. 

Note that, from a field theory point of view it is customary to identify the small 
length scale limit with k — > oo and the large length scale limit with A; — 0. Indeed, 
it is intuitive that there is generically correspondence between s — > and k — > oo or 
between s — > oo and A; — ?• 0. In fact, if one is interested in calculating only these 
limits, it is much easier to suitably truncate the dispersion relation for small or large 
k and then directly calculate the spectral dimension in these limits (as was done in 
reference [10]). 

It is worth mentioning that the correspondence between the two limits (in s and 
k) can be readily seen from the end result of the three examples we gave above. Notice 
that in all three cases the spectral dimension is actually a function of sK n , where K 
is the energy scale suppressing the higher order momentum terms in the dispersion 
relation (and n is the exponent of the lower order term based on dimensional analysis). 
Hence, clearly the limit s — > (s — > oo) will give no different results than the limit 
K —7- (K — > oo). However, given the form of dispersion relations we are considering, 
driving K to infinity is no different than taking k — > 0, and driving K to zero is not 
different than taking k — > oo. 

The general pattern relating the small/large length scales (ultraviolet /infrared lim- 
its) of the spectral dimension suggested by our analytic examples will be be investigated 
in the next two sections using both asymptotic techniques and saddle point methods. 
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3.2 Asymptotic expansions 

Since obtaining explicit analytic results for the spectral dimension seems to be limited 
to very simple dispersion relations, it is worth developing perturbative expansions that 
could help us determine the former from the latter to some desired accuracy in a 
regime of interest. Specifically, if one is interested in studying the very small length 
scale (ultraviolet) or the very large length scale (infrared) behaviour of the spectral 
dimension, then asymptotic techniques are particularly useful. These are based on the 
idea of dividing the exponential appearing in the partition function, e~ s ^ k2 \ into a 
dominant piece plus perturbative corrections. Exactly how one does the division will 
depend on which of the two asymptotic regimes one is interested in probing. 
Consider a dispersion relation of the form 

co 2 = f(k 2 ) = c w k 2w + c w+l k 2w+2 + ■■■ + c z ^k 2z ~ 2 + c z k 2z , (3.23) 

with 1 < w < z. Then one can certainly write both 

exp[-sf(k 2 )} = exp{-sc w k 2w ) {l - sc w+l k 2w+2 + ...}, (3.24) 

and 

exp{-sf(k 2 )} = exp(-sc z k 2z ) {l - sc^k 22 ' 2 + ...}, (3.25) 

by performing expansions around k — and s = respectively. The series in braces 
consist of various positive powers of s and k and are guaranteed convergent for all 
values of s and k. The sub-leading terms appearing in equation (3.24) are all of the 
form s p k 2q with q > p(w + l). Similarly sub-leading terms appearing in equation (3.25) 
are all of the form s p k 2q with q < p(z — 1). 

These two expansions can be very useful for studying the infrared and ultraviolet 
behaviour of the spectral dimension. Moreover, they are not actually limited to poly- 
nomial dispersion relations. Consider instead a general dispersion relation u 2 = f{k 2 ). 
If /(■) does not have a pole at k — 0, then one can always write 

uj 2 = f(k 2 ) = c w k 2w + c w+1 k 2w+2 + ..., (3.26) 

which is sufficient to obtain equation (3.24). On the other hand, as long as /(•) has 
a pole of finite order at k — > oo, one can always write the dispersion relation for very 
large k as a Taylor series of descending powers of k, i.e. 

co 2 = f{k 2 ) = c z k 2z + c z _ x k 2z - 2 + .... (3.27) 

Again, this is sufficient for obtaining equation (3.25). Both requirements are fulfilled 
by reasonable and phenomenologically interesting dispersion relations. 
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We are now ready to proceed to suitable ultraviolet and infrared perturbative 
expansions. Before doing so, it is worth mentioning the following. We have restricted 
our discussion so far to dispersion relations that respect parity invariance. This is a 
very well motivated simplifying assumption: The lowest order term that violates parity 
{uj 2 oc k) is important at low energies and would, therefore, lead to severe disagreement 
with observations. However, one could relax parity invariance and introduce some other 
symmetry that would prohibit the presence of lower-order but not higher-order parity 
violating terms (one example seems to be Hofava's detailed balance condition [18]). In 
any case, even though we will retain parity invariance for the sake of brevity, it is not 
crucial for our approach. After deriving the asymptotic expansions, we will comment 
on how they would be modified if parity invariance were to be abandoned. 

3.2.1 Deep ultraviolet (generalized Seeley— DeWitt expansion) 

We now argue that there is an appropriate generalization of the usual Seeley-De Witt 
expansion beyond ordinary second-order differential operators. If /(—V 2 ) is polynomial 
and contains terms up to order V 2z , then there will exist an asymptotic expansion (in 
terms of positive fractional powers of s) of the form 

ZM = ^ {X>"" + O (^ l «>)}, (3.28) 

ln=0 J 

where we normalize to a,Q = 1, (and C is for our current purposes uninteresting). The 
usual Seeley-De Witt expansion corresponds to z = 1. Note that we have already 
peeled off physical time and are only dealing with the spatial part of the heat kernel. 
To convince yourself that this expansion holds, proceed as follows: Focus on the highest- 
derivative piece. A straightforward computation (using the definition if the T function) 
yields 




(3.29) 



Note that the sub-leading terms are all of the form s p k 2q and so will contribute 

J k d ~ 1 {s p k 2q } exp(-sk 2z ) dk oc s~ dl(2z) ■ s p ■ s~ q/z . (3.30) 

So in particular, for the first sub-leading term in equation (3.25) which is proportional 

to sk 2z ~ 2 , we have 

J k d ~ l {sk 2z - 2 } exp(-sfc 2z ) dk oc s~ d ^ -s-s- 1 - s +l ' z = 8 ~ d ^ ■ s +l ' z . (3.31) 
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Then, at the very least Z(s) will contain the terms 

Z(s) = C{a s- d / {2z ^ + a lS - d ^ ■ s 1/z + ...}. (3.32) 

Higher order terms s p k 2q contribute to the partition function with relative strength 
s p . s -<il z ; that is s^ pz ~ qS} / z . But pz — q > p is a positive integer, so these are all positive 
powers s n / z of s 1 ^. That is, this series for Z(s) (which can at best be an asymptotic 
series) continues with increasing powers of (s 1 ^)™. This implies that Z(s) is indeed 
given by equation (3.28) as claimed. 

Now, turning to the spectral dimension, one can straightforwardly derive 

d s (s) = 1 + - - ^-s l ' z + 0{s 2 ' z ). (3.33) 

So, in the deep ultraviolet (s — > 0) the spectral dimension flows generically (modulo 
our assumptions about parity invariance) to 

d s (s) -»• 1 + (3.34) 

with a prescribed rate of 0(s 1 ^ z ), where 2z is the order of the pole of the dispersion 
relations as k — > oo (and d the topological spatial dimension). We are (for current 
purposes) not particularly interested in the specific value of the coefficient a\, though 
it can be calculated without any particular difficulty When we normalize to a = 1 a 
brief calculation yields 

This quantity a\ is not, in any sense, universal like the usual Seeley-De Witt coefficients 
- ai will depend explicitly on the coefficients in the dispersion relation f(k 2 ). Had we 
allowed for parity violations the only qualitative modification in the spectral dimension 
as given by equation (3.33) would be that z would be allowed to take half-integer values. 

3.2.2 Deep infrared (generalized Feynman expansion) 

Turning to the infrared, we now argue that there is an appropriate "inverted" gen- 
eralization of the Seeley-De Witt expansion beyond ordinary second-order differential 
operators. (This is very close in spirit to the usual Feynman diagram expansion around 
a Gaussian integral.) Specifically, if w = 1 so that the low-momentum dispersion rela- 
tion starts off as f(k 2 ) = k 2 + we assert the existence of an asymptotic expansion 
(now in inverse integer powers of s) of the form 

^) = ^{E^""+ ( s " (iV+1) )V ( 3 - 36 ) 

I n=0 J 
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where we can normalize to do = 1. Even if we do not have Lorentz invariance in the 
low energy limit (w > 1), we can still make the general assertion 

Z & = -pk lt^- n/W + ° (s^ N+1)/w ) ) • (3.37) 

I n=0 J 



To convince yourself that this is true, proceed as follows (many of the technical steps 
are very similar to those for the ultraviolet limit and some details will be suppressed). 
Focus on the lowest-derivative piece, for which we have 



Z ( _ v2)TO (s) = J k d - l e~ sk2w dk oc s~ d/2w . (3.38) 



Furthermore, note that all the sub-leading terms are of the form s p k 2q , and so by a 
minor variant of the previous argument they contribute quantities of relative magnitude 
s p . s -q/w £ the partition function Z(s). But the first sub-leading term in equation 
(3.24) is proportional to sk 2w+2 , so at the very least Z(s) will contain the terms 

Z(s) = C{a s~ d/{2w) + s ■ a lS ~ d/{2w) ■ + ...}. (3.39) 

Higher order terms contribute with relative strength s p ■ s~ q l w = s <yPW ~ qS} / w ', but 
from equation (3.24) we know q > p(w + 1) so pq — q < —p is a negative integer. 
So higher order terms contribute with relative strength s~ n ' w , that is, with positive 
integer powers of s~ l / w . Thus, Z(s) is indeed given by equation (3.37) as claimed. The 
corresponding spectral dimension is then 

d s ( s ) = 1 + - + 2ti lS - 1/w + 0{s- 2 l w ). (3.40) 
w 

Thus deep in the infrared (s — > oo), as long as the low-momentum dispersion relation 
is of the form f(k 2 ) = k 2 + k A + (that is w — 1), the spectral dimension will flow to 
d+1, the topological dimension of spacetime, at a rate 0(1/ s). More generally, it will 
flow to 1 + d/w at a rate of 0(s~ 1 l' w ). As before, had we given up parity invariance, w 
would be allowed to take half-integer values. Finally, di is not in any sense universal 
and can be calculated without difficulty Normalizing to So = 1, a brief calculation 
yields 

_ c w+1 r(jg + 1) , s 



(c) r( 



2w> 



- 19 - 



3.2.3 Example: Rational polynomial dispersion relation 

As previously mentioned, non-projectable Hofava gravity models lead to dispersion 
relations that are rational ratios of polynomials in k 2 [24, 25] so that 

< 3 - 42 > 

Since p(k 2 ) and q(k 2 ) are both polynomials in k 2 we can argue that 

-> c z k 2z + 0(k 2z ~ 2 ) as k^oo. (3.43) 

q(k z ) 

Then we have 

f(k 2 ) = c z k 2z + 0(k 2z ~ 2 ) as k -)• oo . (3.44) 
Therefore for the spectral dimension we have 

d s = l + - + 0(s l ' z ). (3.45) 

z 

A similar procedure works in the infrared. At low k any rational ratio p(k 2 )/q(k 2 ) can 
be expanded as a formally infinite power series in k 2 . Ignoring issues of convergence, 
we can then still apply the generalized Feynman expansion to argue for behaviour of 
the form 

d s = l + - + 0{s- 1 ), (3.46) 
w 

where we define w by p(k 2 )/q(k 2 ) — > k 2w + 0(k 2w+2 ) as k — > 0. 
3.3 Saddle point techniques 

We now have good control of both asymptotic limits, which (where they overlap) are 
in complete agreement with the exact results obtained previously. However, in the 
intermediate regime neither asymptotic expansion need be reliable — neither general- 
ized Seeley-De Witt nor generalized Feynman expansions are guaranteed to be useful. 
Fortunately, in this situation saddle point techniques can provide significant and useful 
approximate information. 

Consider any integral of the form 

/+oo 
exp[h(x)\ dx. (3.47) 
-oo 

The saddle point approximation consists of first locating the maximum of the argument 
h{x) by solving 

h'M = 0, (3.48) 
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and then approximating 
J ft 



exp 



1 

h( x *) + 2^"( x *)( x ~ x *y 



dx 



/+oo 
exp 
-oo 

exp [h(x*)} x 



dx 



2tt 



h"(x*)' 



(3.49) 
(3.50) 

(3.51) 



For most purposes (especially when looking at ratios of integrals) it is sufficient to use 
the even more brutal approximation 



exp[h(x)} dx ~ exp [/i(x*)] 



(3.52) 



See for example reference [27], where similar considerations are applied to the stationary 
phase approximation. 

3.3.1 General framework 

Taking the change of variable u = In k we can express the partition function as 



Z(s) = / exp{ud - s f(e zu )} du 



Note that 



{ud - s f{e 2u )}' = d-2s f{e 2u )e 2u , 
and so the location of the saddle point is at 

d = 2sf(e 2u )e 2u = 2sf{kl)kl 



(3.53) 



(3.54) 



(3.55) 



which implicitly defines k*(s). We can rewrite this equation in terms of the phase and 
group velocities as 



d — 2s -7—— k^ — 2s fphasc(^*) ^'group(^*) k^, 



(3.56) 



implying that the saddle point estimate of the diffusion time can be specified as a 
function of the wavenumber 



s(fc*) 



d 



2fpi lase (/Cifc) ^group 



(3.57) 
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We now see 

Z(s) « ex V [{ud - s f(e 2u )}*} = exp{-s f(h*)}. (3.58) 

Therefore, in view of the saddle point condition in equation (3.55), the leading term in 
the saddle point approximation yields for the spectral dimension 

ds(s)*l + 2sf(kl) = l + jl^. (3.59) 

That is, expressing s in terms of fc*, 

d s (h) = l + d Vphas j k *\ +... (3.60) 

^group \ ) 

s{K) = 2vu (ku — (kY¥- (3 - 61) 

* '-'phase V""*/ ty groupV A '*/ rv * 

If desired, the curve ds(s) can be explicitly constructed by parametrically plotting the 
pair {d s {h), s(fc*)}. 

To once more verify our asymptotic results, let us now consider the situation 

fi(Jfe) ~ k w as fc -> 0; and fi(Jfc) ~ fc* as fc — >■ oo. (3.62) 

We then have 

s (fc*) ~ ^ 00 as fc* — 0; and s(fc*) ~ — — > as fc* — )■ oo , 

2wk% w 2zk% z 

(3.63) 

as expected. Additionally 

d d 
ds{k*) ~ 1 H as fc* — )■ 0; and d s {k*) « lH — as fc* — >• oo , (3.64) 

u> z 

but now there is no requirement that and z be integers. So, the infrared and ultravio- 
let limits for ds(s) that we first encountered in several exact examples, and then verified 
for general polynomial dispersion relations using asymptotic expansions, are now seen 
to have even more generality within the framework of the saddle point approximation. 

3.3.2 Example: Polynomial dispersion relation 

Suppose we have a polynomial dispersion relation 

z 

Q(k) 2 = J2 c ak 2a - (3.65) 
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Then 

W*) = v phase(k) = V/Ea , Cafc2Q ; (3.66) 

so the saddle point approximation implies 

V r k 2a d 

d s {h) *l + d fl * ; s(K) = -= — . (3.67) 

La a °aki 2 ac a kf 

Then: 

d s (s -»■ oo) = 4 0)«1 + -; ds(s -»• 0) = d<j(fc* ->• oo) « 1 + -; (3.68) 

which agrees with our explicit limits via the asymptotic analyses — but now also gives 
us approximate information at intermediate regimes. 

3.3.3 Example: Two-term dispersion relation 

It is worthwhile to explicitly consider the simple two-term dispersion relation 

n(k) 2 = c w k 2w + c z k 2z , (3.69) 
since then we have particularly simple formulae 

d(h)^l + d w * + z * ; s(K) = - — ? — . (3.70) 

3.3.4 Example: Bogoliubov dispersion relation 

Consider the specific case of the Bogoliubov dispersion relation written in the form 



m 2 = k 2 i + — , . (3.7i) 



Then the saddle point approximation yields 

In this particular case we can explicitly invert s(fc*) to find fc*(s): 



«W = *"lVsi + 1 - 11 - (373) 
Thence 

, 5(g) . 1 + ^+f(^H. (3.74) 
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That is 

' 2l + + d/fs/f!) 1 

This compact and explicit (though approximate) formula has the correct limits. (ds{s — > 
oo) = while d s (s — > 0) = l + d/2.) This result can in 2 + 1 dimensions be explicitly 
compared to the exact analytic result obtained via the Gaussian error function. 

3.3.5 Example: Surface waves 

Take a two-dimensional liquid-gas interface, (topological dimension d = 2), and consider 
the surface waves. For example, finite-depth ocean waves ignoring surface tension. 
Then the dispersion relation is well known to be 

tt(k) 2 = gkttmh(kh), (3.76) 

where h is the depth of the ocean. Since this an intrinsically non-relativistic system, as 
we have previously argued it is best to simply consider the spectral dimension of space 
(rather than spacetime). The saddle point approximation leads to 

4sinh(2A;*/i) 4 

«S,space(K*) « „ ; „ Wo ,„ M , ,„ I = — 2kji ' (3.77) 



smh(2k*h) + 2k*h 1 



sinh(2fc*h) 



while 



= ; : — ; o7 — \T- ( 3 - 78 ) 

gk*(ta.nh.(kh) + k^hsech (kh)) 



Then in the asymptotic regimes we have 

ds,s P acc(.s ->• oo) = d s (K -»■ 0) w 2; rf 5iSpacc (s ->• 0) = d 5 (A; + ->• oo) w 4. (3.79) 

So in the infrared one recovers the topological dimension, as expected — but in the 
ultraviolet something slightly unusual happens in that the spectral dimension increases. 
This is ultimately due to the fact that for large k we have Jl(fc) oc k 1 ^ 2 , with an exponent 
that is lower than the Q(k) oc k behaviour characteristic of the infrared. 

Equivalently we could work with fixed /c* and consider the effect of varying h. We 
have 

ds, spa ,cc(h -> 0) w 2; d SiSpSLCC (h ->■ oo) w 4. (3.80) 

This implies that for shallow water surface waves the spectral dimension is 2, whereas 
for deep water surface waves the spectral dimension is 4. 

If you include surface tension the dispersion relation is modified and eventually 
at very short distances a different asymptotic behaviour dominates: f2 oc k 3 ^ 2 at very 
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large k. So in this region the surface tension dominated deep water waves exhibit a 
spectral dimension: 



4 Conclusions 

We have analyzed the concept of the spectral dimension in the continuum and focused 
mainly on scales where curvature effects can be neglected. We have argued that in this 
regime the spectral dimension characterizes a differential operator that is used to define 
the (fictitious) diffusion process, or better yet the dispersion relation associated with 
this operator. We exhibited a natural way of assigning a spectral dimension to any 
dispersion relation, and shown how one can in principle invert the relationship. These 
results establish that the spectral dimension acts as probe of the kinematics of a given 
degree of freedom in this regime. 

We have considered some simple examples of dispersion relations for which one can 
analytically determine the spectral dimension. However, for more general (and realis- 
tic) dispersion relations this analytic determination is rarely feasible. To address this 
issue we have developed approximate techniques that allow one to determine the be- 
haviour of the spectral dimension to a desired accuracy in various regimes of interests. 
In particular, we have presented asymptotics expansions (that resemble Seeley-DeWitt 
and Feynman expansions) which allow one to obtain the infrared and ultraviolet be- 
havior of the spectral dimension (in flat space), and a technique to calculate the flow 
of the spectral dimension at intermediate scales using a saddle point approximation. 

As a major application of the theoretical and technical developments presented 
here, we consider the use of the concept of the spectral dimension as a link between 
discrete and continuum quantum gravity theories. Given that in the former the spectral 
dimension is one of the few known observables, defining it and providing the techniques 
for calculating in the latter is an important step in understanding the continuum limit 
of discrete models. In a companion article [15] we shall focus on this issue in more 
detail. 
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